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Abstract 

We study the minimal speed of propagating fronts of convection reaction diffusion equations of 
the form ut + /-i(ft(u)u x = u xx + f(u) for positive reaction terms with /'(0 > 0. The function 4>(u) 
J> ' is continuous and vanishes at u = 0. A variational principle for the minimal speed of the waves is 

o\ ' 

constructed from which upper and lower bounds are obtained. This permits the a priori assesment 

o. 

of the effect of the convective term on the minimal speed of the traveling fronts. If the convective 
term is not strong enough, it produces no effect on the minimal speed of the fronts. We show that 
if f"(u) I \J /'(0) + pL(j)'(u) < 0, then the minimal speed is given by the linear value 2-^/ /'(0), and 
the convective term has no effect on the minimal speed. The results are illustrated by applying 
. ,— i ', them to the exactly solvable case ut + fiuu x = u xx + u{l — u). Results are also given for the density 

5-h ' dependent diffusion case ut + [icf)(u)u x = (D{u)u x ) x + f{u). 

a | 

PACS numbers: 87.23.Cc, 04.20.Fy, 05.60.Cd, 02.60.Lj 
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I. INTRODUCTION 



The reaction diffusion equation u t = u xx + f(u) has been employed as a simple model 

of phenomena in different areas, population growth, chemical reactions, flame propagation 

and others. In the classical Fisher case f (u) = u(l — u) , a front propagating with speed 

c^-pp = 2 joins the two equilibrium points |2j. The time evolution for general reaction terms 

was solved by Aronson and Weinberger [AW] who showed that sufficiently localized initial 

conditions evolve into a front which propagates with speed c* such that 2^J f'(0) < c* < 

2^/sup(/(w) /u). The asymptotic speed of propagation is the minimal speed for which a 

monotonic front joining the stable to unstable equilibrium point exists. Existence proofs 

give limited quantitative information on the dependence of the speed of the front on the 

parameters of the problem For this reason different variational methods have been 

developed. For the one dimensional case it has been shown that this minimal speed can be 

rn 

derived either from a local variational principle of the minimax type [a], or from an integral 
variational principle Minimax variational principles for the speed of fronts in several 

dimensions and for inhomogeneous environments have also been established jj, ^. 

In many processes, in addition to diffusion, motion can also be due to advection or convec- 
tion. Nonlinear advection terms arise naturally, for example, in the motion of chemotactic 
cells. In a simple one dimensional model, denoting by p the density of bacteria, chemotactic 
to a single chemical element of concentration s(x, t) the density evolves according to 

p t = [Dp x - pxs x ]x + f(p), 

where diffusion, chemotaxis and growth have been considered. There is some evidence 
that, in certain cases, the rate of chemical consumption is due mainly to the ability of the 
bacteria to consume it. In that case 

s t = -kp, 

where diffusion of the chemical has been neglected (arguments to justify this approximation, 
together with the choice of constant D and x are given in j^). If we now look for traveling 
wave solutions s = s(x — ct), p = p(x — ct), then s t = —cs x , therefore s x = kp/c, and the 
problem reduces to a single differential equation for p, namely, 

Pt = D Pxx -^(pX + f(p). (I) 
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The more elaborate models of Keller and Segel for chemotaxis Q], which include diffusion 
the chemical and other effects, have been considered to explain chemotactic collapse 
ll^ and references therein) and other phenomena. The derivation of these equations 
from transport theory and the assumptions involved in them have been studied recently [13] . 
In addition to these biological processes, equations analogous to (JI} appear when modeling 
the Gunn effect in semiconductors and in other physical phenomena 3> (3 • Equation (JTJ) 
for a Fisher type reaction term f(u) = u(l — u) has been studied in [16]. An extensive study 
of the existence of traveling waves of nonlinear diffusion-reaction-convection equations which 



mc 
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udes a review of many results to which we refer for additional references is contained in 



In this work we concentrate on the equation with a general convective term which, suitably 
scaled, we write as 

ut + n<p(u)u x = u xx + /(«), (2) 

where the reaction term f(u) is a continuous function with continuous derivative in [0, 1] 
and satisfies 

/(0) = /(l) = 0, /'(0)>0, /'(1)<0 and />0 in (0,1). 

The function <p{u) is a continuous function with continuous derivative in [0, 1]. Without loss 
of generality we may assume that in addition 0(0) = 0, since otherwise only a uniform shift 
in the speed is introduced. The parameter \i is positive. 

For equation (2), the existence of monotonic decaying traveling fronts u(x — ct) for anj 
wave speed greater than a critical value c* has been proved recently [18j . Moreover, in 
the following estimate for the threshold value c* is obtained, 



2Jf'(0) < c* < 2 J sup + max ud)(u). (3) 

« e( o,i] u «e[o,ir 



Analogous results for density dependent diffusion are also established in The conver- 



gence of some initial conditions to a monotonic traveling front has been proved for 
systems in which the minimal speed is strictly greater than the linear value ci = 2a//'(0). 

We will show that the minimal speed c* for the existence of a monotonic decaying front 
u(x — ct) joining the stable equilibrium u — 1 to the unstable equilibrium u = obeys the 
variational principle 

c* = sup£(g), (4) 



with 

£{g) = J 1 (2 y /f(u)g(u)(-g'(u)) + fuf>(u)g(u)) du 

jlg{u)du 

and the supremum is taken over the set S of all positive, monotonic decreasing functions 
g(u) for which the integrals in (jSJ) exist and g(l) = 0. From here it will follow that 

2V? 7 (0) = c L <c,< inf sup (a + -f{u) + fM<j>(u)). (6) 

a>0 «e[o,i] a 

From the variational expression (j3J) one may obtain the value of the minimal speed with any 
desired accuracy, and the inequalities (jBJ) enable us to characterize the reaction terms for 
which the speed is the linear value cl- More precisely, if 

f"( u ) , 

J-^L + ^'iu) <0, then c = 2v7W 

V / (°) 

The bound (JHJ) is also derived from the variational principle. The generalization to density 
dependent diffusion is given as a direct extension of the previous results. 



II. MINIMAL SPEED OF TRAVELING FRONTS 

Traveling monotonic decaying fronts u(x — ct) of (0) satisfy the ordinary differential 
equation 

u zz + (c — /i(j)(u))u z + f(u) = lim u — 1, lim it = 0, u z < 0, 

where z = x — ct. It is convenient to work in phase space; defining as usual p(u) = —u z , the 
problem reduces to finding the solutions of 

P(«) ^ " (c - A*^(«))p(u) + /(«) = 0, (7) 

with 

p(0)=p(l) = 0, and p > 0. 

We first perform the linear analysis around the endpoints u = and u = 1 which may 
provide restrictions on the allowable speed. These results will also be needed when proving 
the existence of a variational principle. Near u — 0, p(u) = mu where m is the larger root 
of m 2 — F(0)m + /'(0) = 0. For convenience, we have defined F{u) = c — fJ,<j>(u). This root 
is given by 

m 



F(0) + " 4/'(0) 
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The condition that m be real imposes the restriction F 2 (0) > 4/'(0). Written explicitly this 
bound is, 

c>2^f{0) = c L . (8) 



Near u = 1, p = r(l — u) where r is the positive root of r 2 + F(l)r + f'(l) = 0, namely, 

_ -F(l) + v/F2(l)-4f(l) 



(9) 



2 

No additional restriction on the range of allowable speeds is imposed from the expression 
above, since by hypothesis /'(l) < 0. 

In addition to the linear constraint (JBJ), a simple constraint is found from direct integration 

of Eq.(JZJ). Dividing by p(u) and integrating between and 1, we have 

pi pi 

c=u 6(u)du + / . du. 
Jo Jo PW 



Since / and p are positive in (0, 1), we obtain 



Out 4>{u)du. 
Jo 

A. Variational principle 

In this section we construct a variational principle from which the exact speed of the 
front may be calculated. Let g be any positive function in (0,1) such that h = —dg/du > 0. 
Multiplying Eq.(J7J)) by g/p and integrating with respect to u we find that 

c [ g du = [ (hp+— g)du + a [ (pgdu 
Jo Jo V P J Jo 

where the first term on the right-hand side is obtained after integration by parts. However 
since p, h, f, andg are positive, we have that for every fixed u 

fg 

np-\- 

so that, 



p 



c > fiPVTgK + rigVu = £{g) (lQ) 

Jo 9du 

To show that this is a variational principle we must prove that there exists a function g = g 
for which equality holds. Equality is attained for g = g such that 

i gf 

ph = -pg = — . 

P 
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Using Eq.(J7j) to eliminate f(u) we have that 

g'{u) p'{u) F{u) 



g(u) p(u) p(u) 

which can be integrated to obtain 



g(u) = p(u) exp 



with < m < 1- (11) 



Since p vanishes at and 1, we must analyze the behavior of g at these points in order to 
ensure the convergence of the integrals in fjlUjl . 

At u — 1, since p(u) = r(l — u), and since F{u) is continuous at 1, we obtain that 
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r (l_ u )(i+m)/r). 



From the expression for r, Q we see, that, since /'(I) < 0, for any value of F(l) the 
exponent 1 + F(l)/r is positive, hence g(l) = 0. 

Near u = 0, since p = ma and since F is continuous at zero, we find that 

m 
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u (F(0)/m)-l ' 

The integrals in (fTU|) converge provided 

F(0) 

- 1 < 1. 

m 

This condition is satisfied whenever C* > c^. That is, whenever c* > there exists a 
function <? = g for which equality holds in (fTUj) or , equivalently, c* = max £((7) = £(g). 

On the other hand, when c* = cl, there does not exist a function g in the set 5 of 
admissible functions for which equality holds in (jlOJl . Consider however the trial functions 

c/ a (u) = m"- 1 - 1. 

Clearly g a E S for any < a < 1. Moreover one can check that (see the Appendix) 

\im£(g a ) = c L = c*, 



therefore 



c* = sup^(^) 
36-s 



in this case. Notice that in this last case the maximum is not attained since the limiting 
function go does not belong to S. This concludes the proof of our variational principle. 
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B. Upper and Lower Bounds 



The variational principle provides lower bounds with suitably chosen trial functions, which 
can be arbitrarily close to the exact value of the speed. The fact that it is a variational 
principle for which equality holds, enables one to obtain also an upper bound to the speed. 

To obtain an upper bound we use the fact that 

2ab < aa 2 + -b 2 , 

a 

Then the following inequality holds: 



with a > 0. 



(12) 



I^Jgh = 2g^]hfg <g(a + 

V a 9 J 



where we used the inequality above with a = 1,6= y/ fh/g. 
Then, 



c* = sup 

9 



lo^Vf 9 h + y<t>g)du 

Jo 9 du 



< sup 

9 



Jo g[& + fh/(ag) + n<j>] du 
Jo 9 du 



The second term in the last expression can be integrated by parts. The boundary term fg\\ 
vanishes and we obtain 

Jo g(a + f'/a + u<j))du 



< sup 

9 

< sup 

«e[o,i] 



Jo9du 



a + 



a 



(13) 



The above inequality holds for any positive a, hence 



c* < inf sup 

a uG[0,l] 



(14) 



a H — f'Ou) + u(p(u) 
a 

The bound (|14p in the case u = differs from the classical Aronson- Weinberger Q result 
for fronts of the parabolic reaction diffusion equation c < caw = 2sup Mg [ 01 ] vTR/^. The 
bound © obtained in [18] on the other hand, reduces, when /x = to the classical Aronson- 
Weinberger result. Here we show that this last bound can be derived from the variational 
principle as well. Using the inequality (JT2J), now with a = \J fg/u and b = y/uh, we have 
that 

fg 



2 I \fjghdu < 
'o 



a 1 — nu 

u a 



du 



1 

a— + -g 

u a 



du, 



where the last expression is obtained after integrating the second term by parts. We have 
then 

c* < sup — i < sup 

g J gdu «€[o,i] 

Choosing a = 1/ sup a/77" we obtain 

c* < sup 
«e]o,i] 

In the classical AW case /i = 0, we know that when the reaction term is concave then 
sup2\ v //(u)/u = 2a//'(0). In this case the upper bound coincides with the linear lower 
bound Cl and the minimal speed is univocally determined. A similar criterion can be ob- 
tained in the present problem. The minimal speed for the existence of a front is known 
unambiguosly to be the linear value whenever the upper bound ()14j) coincides with the lin- 
ear lower bound cl- A sufficient condition for this to occur is that the supremum of the 
function K{u) = a + f'(u)/a + n<$>(u) in (0, 1] does not exceed the value of K at the origin. 
Effectively, if sup n A'(w) = K(0) = a + /'(0)/a; minimizing with respect to a we obtain 
a = y/f'(0) and the upper bound is precisely the linear value. A sufficient condition that 
guarantees that the maximum (supremum) of K occurs at zero is that K(u) be decreasing. 
With a = \J /'(0) this condition is 

^^L+/i0'( M )<o. 

VfW) 

Whenever this condition is fulfilled, for all u, we know that the minimal speed of a monotonic 
front is the linear value cl = 2a//'(0). Again as it occurs in the standard case fi = 0, this 
condition is sufficient but not necessary. 



a — I h \i<p 

u a 



u 



+ u max <b(u) 
«e[o,i] 



III. AN EXACTLY SOLVABLE CASE 



' iere we illustrate the above results by applying them to the exactly solvable case discussed 
in a Fisher type reaction term f(u) = u(l — u) and the simplest convective term 

4>(u) = u. By means of a phase space analysis Murray found that the minimal speed 



for the existence of a monotonic decaying front is 

//, > 2 

(15) 



+ f if /x > 2 



if n < 2 



Here we show that the results of the previous section allow the exact determination of the 
speed. 

In this example the linear marginal stability value is given by cl = 2. We first use the 
variational principle to obtain a lower bound. Take the trial function 

g(u) = ^ - with < A < 1. 

A straightforward integration of (|10p leads to 

c>2VA + |(l-A) =c(A). 

If fi > 2 the maximum of c(A) occurs for A = 4//i 2 and it is given by 2//i + /x/2. For /i < 2, 
however, the supremum of c(A) occurs as A — > 1. We have then 

c* > sup c(A) = - + for /i > 2, 

/! 2 

and 

c* > sup c(A) = 2 for /x < 2. 

To obtain an upper bound we use (fT3j) . that is, 

1 2 

c* < sup (cH \- u(fi )) Vet > 0. 

ue[o,i] a a 

We will separate the two cases, fi < 2 and /i > 2. 
If /i < 2, choose a = 1, then 

c* < sup (2 + «(//- 2)) = 2. 

«6[0,1] 

If > 2 choose a = 2/fi, then 

c* < sup - + £ = - + 

«6[0,1] V/ 1 ^/ ^ Z 

The lower bound obtained from the variational expression coincides with the upper bound 
obtained from (jUty , therefore we know with certainty that the minimal speed is indeed, (II 5j) 
which had been previously demonstrated by phase space methods. 

Notice that constitutes a poorer bound in this case. Effectively, from Q it follows 
that c* < 2 + yU. 
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IV. DENSITY DEPENDENT DIFFUSION 



The effect of the convective term on the minimal speed of fronts of the reaction diffu- 
sion equation for non constant difffusion follows in a simple way from the previous results. 
Consider traveling fronts of the equation 

u t + n <p(u) u x = (D(u)u x ) x + f(u), 

where f(u) and 4>(u) satisfy the properties spelled in the previous sections. The diffusion 
coefficient D{u) is continuous and Diu) > in (0,1]. -D(O) is either positive or zero. By 



a suitable change of variables [20, |2l| the equation for the fronts is reduced of the usual 
reaction diffusion equation with a reaction term f(u) = D(u)f(u). This reaction term 
satisfies / > 0, and /'(0) = _D(0)/'(0). We must distinguish two cases. If D(0) ^ 0, then 
/ satisfies the same properties as / and the results of the previous sections can be applied 
directly. If D(0) = 0, then f'(0) = and we expect a sharp wavefront. In this case it has 
been shown J3|l that the front approaches u = as cu. A variational principle exists also 
in this case [22] . We have then that, in both cases, the minimal speed of the wavefronts is 
given by 



sup 

9 



jh2y/D(u)f(u)g(u)(-g'(u))+ ^(u)g(u)) du 



Jo 9(u) du 

where the supremum is taken over all positive monotonic decreasing functions g(u) for which 
the integrals exist and g(l) = 0. Upper and lower bounds can be obtained following the 
methods of the previous sections. We do not spell them out here. 



V. SUMMARY 



We have studied the effect of a convective term on the speed of monotonic reaction- 
diffusion fronts. The minimal speed for the existence of fronts has been showed to derive 
from a variational principle from which the exact speed can be determined in principle. 
The existence of this variational characterization permits the obtention of upper and lower 
bounds. The classical result that establishes that for concave reaction terms, the minimal 
speed of the fronts is the linear or KPP value is extended to the case where convective terms 
are present. The extension to the the case of density dependent diffusion has been given for 
positive diffusion terms. 

10 



We have found that a convective term increases the minimal speed of the traveling front 
only if it is sufficiently strong, if not, the minimal speed is determined by the reaction term 
alone. 
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APPENDIX A 

In this appendix we prove that 

lim^(^)=c L = 2 A /7 7 (0), 

a^0 

where g a = u a ~ x — 1 with < a < 1. 

Since g a (u)du = (1 — a) /a, and h a — (1 — a)u a ~ 2 , we may write 

£(g a ) = Ji(a) + J 2 (a), 

where 

Jl(«) = ~^L= [ y/f( U )(u*°<-3- U «-2)du, 
Vl - OL JO 

and 

JJ a ) = / <t)(u)(u a - 1 - l)du. 

Since 0(0) = and <p is continuous the integral in J 2 has a finite value when a = 0. Then, 
due to the overall multiplicative factor of a, we see that 

lim J2(a) = 0. 

To show that lim J\(a) = 2a//'(0), as a — > 0, write 



2a 



Ji(a) = ^^/" vV/'tOXw 2 "- 3 - u a ~ 2 ) du + K(a), 
Vl - a Jo 



where 

2a 



K(a) 



VT 



a 



[ V 'f(u)(lL 2a ~ 3 - U a - 2 )du - [ y/uf (0)(u 2a ~ 3 - U a ~ 2 )du 

Jo Jo 
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The first integral is 



20 : /V ~/'(0)(u^-U-') rfu- ^ 



3+a \ • 



VI - a Jo y/l-a r(2^iy; 

Now we prove that lim Q ^ K{a) = 0. 

Vl - a Jo 
But j-^/a — \/&| < y/\b — a[, therefore 

VI — ct Jo 

Since /(m) and its derivative are continuous, in [0, 1], there exist d > 0, g > such that 

!/(«)- «/'(o)| 



< du q . 



u 



In particular, if f(u) is analytic in a neighborhood of 0, q — 1. Using this inequality in the 
expression above, we have that 

On/ 

\K(a)\ < . / y /dui+ 1 (u 2a - 3 -u a - 2 )du. 
Vl - a Jo 

Finally, since u 2a ~ 3 — u a ~ 2 < u 2a ~ 3 , 



Vl - a Jo 



Therefore, lim a ^o — 0. 

To sum up, 



„-„ ' a-,0 V /T^ r(g^) 
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